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I. INTRODUCTION 


The mechanism of gravothermal instability, discovered by Antonov is an important 
phenomena in gravitational thermodynamics. It has been very helpful for an extensive 
research concerning statistical mechanics of long range interactions systems in several helds 
in physics . This connection with thermodynamics and statistical mechanics has motivated 
us to investigate statistically the gravothermal catastrophe 

At the same time, there are theoretical evidences that the understanding of gravity has 
been greatly benefited from possible connection with thermodynamics. Pioneering works of 
Bekenstein and Hawking [ 3 ] have described this issue. For example, quantities as area and 
mass of black-holes are associated with entropy and temperature respectively. Working on 
this subject, Jacobson jsf interpreted Einstein field equations as a thermodynamic identity. 
Padmanabhan gave an interpretation of gravity as an equipartition theorem. 

Recently, Verlinde l3| brought an heuristic derivation of gravity, both Newtonian and 
relativistic, at least for static spacetime. The equipartition law of energy has also played an 
important role. On the other hand, one can ask: what is the contribution of gravitational 
models in thermodynamics theories? 

The concept introduced by Verlinde is analogous to Jacobson’s j^, who proposed a ther¬ 
modynamic derivation of Einstein’s equations. The result has shown that the gravitation 
law derived by Newton can be interpreted as an entropic force originated by perturbations 
in the information “manifold” caused by the motion of a massive body when it moves away 
from the holographic screen. Verlinde used this idea together with the Unruh result 0 and 
he obtained Newton’s second law. The entropic force together with the holographic prin¬ 
ciple and the equipartition law of energy resulted in Newton’s law of gravitation. Besides, 
Verlinde’s ideas have been used in cosmology [l^ . 

Currently, two of the most investigated extensions of the usual Boltzmann-Gibbs (BG) 
theory are Tsallis thermostatistics theory [l3l-[l5l| (TT) and Kaniadakis power law statistics 
UMi (KS). The former case, TT, initially considers the entropy formula as a nonextensive 
(NE) quantity where there is a parameter q that measures the so-called degree of nonexten- 
sivity. This formalism has been successfully applied in many physical models. An important 
feature is that when g —)■ 1 we recover the usual BG theory. On the other hand, the so called 
KS naturally emerges from the context of special relativity and in the kinetic interaction 
principle. This formalism has also been successfully applied in many physical models. There 
is a parameter k in the KS that in the limit k —)■ 0 the BG theory is also recovered. 

This paper has two parts: in the hrst one, the self-gravitating system jl^ (and references 
therein), is discussed in the context of the TT and KS formalisms. We found that for both 
formalisms the dependence of the specihc heat, Cv^ on the nonextensive parameter q and 
on the K parameter give rise to a negative branch, a result that features the gravothermal 
collapse. In the second part, we have used the TT formalism and the KS in the Verlinde 
formalism 1^. As a result, the equipartition theorem is derived in the framework of KS 


leading to a modihed gravitational constant. 

This paper is organized in as follows: In section [TTl the main steps of Tsallis, Kaniadakis 
and Verlinde’s formalisms are reviewed. In section im a modified gravitational constant is 
obtained in the light of Kaniadakis’ formalism. In Section HVl in an holographic background, 
we have obtained relations between q and k. In section |V]the gravothermal collapse is studied 
in the light of TT and KS formalisms. In section the relative errors of the gravitational 
constants are calculated . The conclusions and perspectives are presented in the last section. 
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II. NONEXTENSIVE STATISTICS AND HOLOGRAPHIC ENTROPY 


The study of entropy has been an interesting task through recent years thanks to the 
fact that it can be understood as a measure of information loss concerning the microscopic 
degrees of freedom of a physical system, when describing it in terms of macroscopic vari¬ 
ables. Appearing in different scenarios, we can conclude that entropy can be considered as 
a consequence of the gravitational framework @,0- These issues motivated some of us to 
consider other alternatives to the standard BG theory in order to work with Verlinde’s ideas 
together with other subiectsfl0. 

The objective of this section is to provide the reader with the main tools that will be 
used in the following sections. Although both formalisms are well known in the literature, 
these brief reviews can emphasize precisely that there is a connection between both ideas 
which was established recently j^ . 


Tsallis’ formalism 


Tsallis [13| has proposed an important extension of the BG statistical theory and curi¬ 


ously, in a technical terminology, this model is also currently referred to as nonextensive 
statistical mechanics. TT formalism defines a nonadditive entropy given by 


Z^i=l fi 


Sq = ks 


q 


w 

E 

2=1 


P, 


= 1 


( 2 . 1 ) 


where Pi is the probability of the system to be in a microstate, W is the total number of 
conhgurations and g, known in the current literature as Tsallis parameter or NE parameter, 
is a real parameter which quantihes the degree of nonextensivity. The definition of entropy 
in TT formalism possesses the usual properties of positivity, equiprobability, concavity and 
irreversibility and motivated the study of multifractals systems. It is important to stress 
that Tsallis formalism contains the BG statistics as a particular case in the limit g —>■ 1 


where the usual additivity of entropy is recovered. Plastino and Lima [2l| have derived a 


NE equipartition law of energy. It has been shown that the kinetic foundations of Tsallis’ 
NE statistics lead to a velocity distribution for free particles given by 


f,iv) = Bq 


l-(l-g) 


mv 


1/1-9 


2,k]^T 

where Bq is a normalization constant. The expectation value of is given by 


( 2 . 2 ) 


< > g = -^0 


/o°° fqV^dv 
2 ksT 


(2.3) 


5 — 3g m 

The equipartition theorem is then obtained by using 
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(2.4) 


1 9 

Eg = -Nm < V >q, 

where we arrive at 

Eg = ^ NksT . (2.5) 

b — 6q 

The range of g is 0 < g < 5/3. For q = 5/3 (critical value) the expression of the equipartition 
law of energy, Eq. fl2.5l) . diverges. It is also easy to observe that for g = 1, the classical 
equipartition theorem for each microscopic degrees of freedom can be recovered. 


B. Kaniadakis’ statistics 


Kaniadakis statistics [16|, also called K-statistics, similarly to TT formalism generalizes 


the standard BG statistics initially by the introduction of ^-exponential and ^-logarithm 
dehned by 


exp^if) 


(yi + K^p + ft/) 


K 

1 


( 2 . 6 ) 


fK, £ — K 

ln«(/) = (2.7) 

with the following operation being satished 

ln«, (exp^if)) = exp^ = /• (2.8) 

By Eqs. (12.bh and (I2.7p we can observe that the ^-parameter deforms the usual dehnitions 
of the exponential and logarithm functions. 

The K-entropy associated with this ^-framework is given by 


SM) 


d^pf 


r - n 

2k 


(2.9) 


which recovers the BG entropy in the limit k —)■ 0. It is important to mention here that the 
K-entropy satished the properties of concavity, additivity and extensivity. Tsallis’ entropy 
satishes the property of concavity and extensivity but not additivity. This property is not 
fundamental, in principle. The ^-statistics has been successfully applied in many experi¬ 
mental fronts. As an example we can mention cosmic rays [1^, quark-gluon plasma j^ . 
kinetic models describing a gas of interacting atoms and photons |25| and hnancial models 

13 . 

The kinetic foundations of ^-statistics lead to a velocity distribution for free particles 
given by 


/k(^) 



1 

Kmv^ 


2kBT 


( 2 . 10 ) 
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The expectation value of is given by 


< >^= 


/o°° Jk v'^dv 

ir 


Using the integral relation 


17 


( 2 . 11 ) 


we have 


r—1 

X expt, 


-x)dx 



( 2 . 12 ) 




The K-equipartition theorem is then obtained by using 


= -Nm < >^, 


( 2 . 13 ) 


( 2 . 14 ) 


where we arrive at 


2 (i-|«))2(i + ifc)2K r(^ + f)r(^-i) 


(2.15) 


The range of«)isO < k < 2/3. For k = 2/3 (critical value) the expression of the equipartition 
law of energy, Eq. fl2.15p . diverges. For k = 0, the classical equipartition theorem for each 
microscopic degrees of freedom can be recovered. 


Verlinde’s Formalism 


The formalism proposed by E. Verlinde [10 


obtains the gravitational acceleration by 
using the holographic principle and the well known equipartition law of energy. His ideas 
relied on the fact that gravitation can be considered universal and independent of the details 
of the spacetime microstructure. Besides, he brought new concepts concerning holography 
since the holographic principle must unify matter, gravity and quantum mechanics. 

The model considers a spherical surface as being the holographic screen, with a particle 
of mass M positioned in its center. The holographic screen can be imagined as a storage 
device for information. The number of bits, which is the smallest unit of information in the 
holographic screen, is assumed to be proportional to the holographic screen area A 


t' D 


( 2 . 16 ) 
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where A = 47rr^ and Ip = \J~^ is the Planck length and /p is the Planck area. In Ver- 
linde’s framework one can suppose that the bits total energy on the screen is given by the 
equipartition law of energy 


E = ]^NkBT. (2.17) 

It is important to notice that the usual equipartition theorem in Eq. fl2.17p . can be derived 
from the usual BG thermostatistics. Let us consider that the energy of the particle inside 
the holographic screen is equally divided by all bits in such a manner that we can have the 
expression 


MG = - NksT. 
2 


(2.18) 


With Eq. fl2.16p and using the Unruh temperature equation llj given by 


, ^ 1 ha 

ksT = - -, 

27r c 


(2.19) 


we are able to obtain the (absolute) gravitational acceleration formula 


llc^ M 

h 

G^. ( 2 . 20 ) 

From Eq. fl2.20p we can see that the Newton constant G is just written in terms of the 
fundamental constants, G = /h. 

III. MODIFIED GRAVITATIONAL CONSTANT 

As an application of NE equipartition theorem in Verlinde’s formalism we can use the NE 
equipartition formula, i.e., Eq. fl2.5p . Hence, we can obtain a modihed acceleration formula 
given by Q 

^ M , , 

where Gq is an effective gravitational constant which is written as 

( 3 . 2 ) 

From result (13.2p we can observe that the effective gravitational constant depends on the NE 
parameter q. For example, when g = 1 we have Gq = G (BG scenario) and for q = 5/3 we 
have the curious and hypothetical result which is Gq = 0. This result shows us that g = 5/3 
is an upper bound limit when we are dealing with the holographic screen. Notice that this 
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approach is different from the one demonstrated in js^, where the authors considered in 
their model that the number of states is proportional to the volume and not to the area of 
the holographic screen. 

On the other hand, if we use the Kaniadakis equipartition theorem, Eq. fl2.15p . in the 
Verlinde formalism, the modified acceleration formula is given by 

a = G, ^ , (3.3) 

where is an effective gravitational constant which is written as 


(1 - i«:)2 (1 + 3,) r(^_3)r(^+i) ’ 

which can be written alternatively as G = f{K)Gi^ used in Eq. (4.10)-(4.13), as we have 
explained in Section IV. 

From result 03.41) we can observe that the effective gravitational constant depends on the 
K parameter. For example, using the Gamma functions properties written in Eq. 05. 8 p we 
obtain for k = 0 that G^ = G (BG scenario). For k = 2/3 and taking into account that 
the only Gamma function that diverge for this value (r(^ — |)) is in the denominator of 
Eq. 03.41) . we have a result which is G^ = 0. 

From the limits k = 0, g = 1, which we recover the BG statistics, and k = 2/3, q = 5/3, 
we can established a linear relation between k and q written as 


K = q — 1. (3.5) 

In particular, the limits k = 0, g = 1 and k = |,g = | lead to Gq = G^. It is clear that 
relation (13. 5 p is valid only for the range 1 < g < 5/3. In Figure 2 we have plotted Gq/G, 
Eq. (13.21) . and G^/G, Eq. (l3.4p . as function of ^-parameter. The dashed line represents the 
Tsallis formalism and the solid line represents the Kaniadakis formalism. 

In order to explicit the difference between the gravitational constants obtained by Tsallis, 
Eq. 03.2p . and Kaniadakis (Eq. 03.4p i formalisms, we define a new function as 

DG = (3.6) 

Cjr 

In Figure 3 we plot DG as a function of the deformation parameter k. 


IV. AN ALGEBRAIC ALTERNATIVE WAY OF OBTAINING THE RELATION¬ 
SHIP BETWEEN q AND k 


We know that to establish values and/or intervals for g is important in different phys¬ 
ical scenarios, as described in the last sections (and in the quoted references). Hence, we 
believe that it would be interesting to obtain g -relations. Having said that let us make an 
assumption and test its validity. In other words, let us explain that the relations obtained 
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FIG. 1; The normalized modified gravitational constants for Tsallis (dashed line= Gg/G) and 
Kaniadakis (solid line= formalisms as functions of K.Note that Gg is closer to G than G^- 

Also note that the strength for gravitational field is greater in TT formalism than in the Kappa 
one. 


from now on will be obtained in the entropic (Verlinde) scenario explored and it is now a 
general result since TT and KS are very different formalisms. So, what we will see now 
can be realized as an algebraic conjecture and not an equivalence between both formalisms. 
These words have to be well understood. 

From fl3.2p and fl3.4l) let us assume that 


2 _ (i-F)Mi + k) r(v-|)r(v + i) 

5-3, (l-|K)Ml+i«:)2K r(i + f)r(i-i)' 

and for k —> 0 we must have q —> 1. Using Eq. fIS.Sp we have that 


2 _ (i-h) (i + f*^) 


(4.2) 


Since k is very small {k —)■ 0), let us eliminate k 3-terms. So, from fl4.2p we have that 


3 5 

L2(5-3g) ^ 4, 


5 1 

L5-3g ^ 2, 


k + 


5-3g 


1 = 0 , 


(4.3) 


and using /c = 0 in f|4.3p we have that q = 1. From fl4.2|) again, if we eliminate /c^ and 
fc^-terms, we have that 
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FIG. 2: The normalized difference between modified gravitational constants obtained by Tsallis 
and Kaniadakis statistics, Ea. (l3.6l) . as a function of k. 


2 5-3gL 2 . 


(4.4) 


so /c = 0 ? = 1- 

From fl4.3p . let us invert the order , i.e., let us make g = 1 in fl4.3p we have that = 

0 All = 0 and k 2 = 3/2. And the root that matters is /c = 0. From (14.2p we have that 


+ - 


2/31 - 15g\ 


y l + 3q l + 3q 


2 4(5 - q) ^ 8(1 - q) 


l + 3g 


= 0 , 


and again, if we make g = 1 in (14.51) we have that 

k{U‘ + - 4 ) =0, 


(4.6) 


(4.6) 


and the roots are 


ki = 0 , 

k2 = V^), 

fcs = -^(2 + x/T3), (4.7) 

and the root fc = 0 is the one. From these results we can see that the assumption in Eq. (14.ip 
is a new relation between k and g (concerning Verlinde’s scenario) that conhrms k = q — 1. 
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With these results we can write the expressions using KS we have to use that G —> fNsG 
where in general terms we have that 


for TT, 
~ \ fin) for KS . 


(4.8) 


In the near future we will use this compact form to relate the gravitational constants. 


V. NEGATIVE HEAT CAPACITY IN SELF-GRAVITATING SYSTEM AND 
GRAVOTHERMAL COLLAPSE 


The gravothermal instability described by Antonov jl[ has an important role in gravi¬ 
tational thermodynamics [2l-[^, namely, the thermodynamics of self-gravitating systems. A 
self-gravitating system is a system where its constituents interact with each other trough 
gravitational forces. 

The discovery of Antonov analyzes a self-gravitating system conhned inside a uniform 
quadratic box. In this scenario, no maximum entropy state can exist below a critical energy. 
After that, this result for energy below this critical value was obtained, where the system 


28| . This is the “gravothermal catastrophe 


” or “Antonov’s 


would collapse and overheat 
instability.” It is related to the very specihc property of self-gravitating systems to have 
negative specihc heats 3^. In other words, when discussed in the context of micro-canonical 
ensemble a self-gravitating system exhibits a negative heat capacity [^. The negativeness 
of the heat capacity has been extensively studied in many astrophysical phenomena. For 
example, we can mention globular cluster and black hole thermodynamics j3l|. 

The gravothermal catastrophe plays a relevant role in the evolution of globular clusters. 
Concerning dense clusters of compact stars, the gravothermal catastrophe can originate 
the formation of supermassive black holes. And statistical mechanics is also relevant for 
collisionless self-gravitating systems 


32 


A. Negative heat capacity 

A negative heat capacity Cy can be explained through the virial theorem for inverse 
square forces concerning bounded systems. It connects the average kinetic energy (K) with 
the average potential energy (V) 


and the total energy is given by 


<K > = 


— < 1 / > 

2 


(5.1) 


E =< K > + <V > (6.2) 

and using fl5.ll) one has that 

E = - <K> . (5.3) 

For moving particles we have that K = ^NKbT and consequently the specihc heat is 
negative 
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= f = (5.4) 

where, for systems with negative Cy that enters in contact with large thermal reservoir, will 
have fluctuations. These fluctuations add energy and make its transient temperature lower. 
It causes inward heat flow which will head it to even lower temperatures. Consequently, 
— Cv systems cannot have thermal equilibrium . 

The self-gravitating system, according to the virial theorem, has the total energy E = —K 
where K is the kinetic energy, K = ^NksT. Therefore, the heat capacity of the system is 




(5.5) 


which is a negative quantity. This negative value indicates an unusual behavior in which the 
system grows hotter. Silva and Alcaniz by using the kinetic foundations of TT, found 
a relation for the heat capacity as 


Cl = (5.6) 

which reduces to the Maxwellian limit Cy = — for q = 1. The negative value of 

Eq. fl5.6p shows that the particular g-dependence of Cy leads to the classical gravothermal 
catastrophe (Cy < 0) for the range 0 < g < 5/3. 

In the K-formalism, k is the equipartition formula in a 3-dimensional space given by Eq. 
fl2.15p and the heat capacity (Cy = ^) in the k- formalism is given by 


(1-HMi + H r(A-|)r(A + i) 
2 {i-lK)Hi + lk)2K r(A + |)r(A-i) 


(5.7) 


Since 


r( 


2k, 


Dr( 


2k 


+ i) 


lim , , , 

‘-»r(i + f)r(i-i) 


= 2k, 


(5.8) 


the Maxwellian limit, Cy = — ^Nks, is determined for k = 0. The normalized heat capac¬ 
ities of Tsallis’ formalism, Cy/Nks, (dashed line), and Kaniadakis’ formalism, Cy/Nks, 
(solid line), as functions of k, are plotted in Figure 1. As we can see, the heat capacity in 
the K, formalism is negative in the range 0 < k < 2/3. Thus, similarly to TT formalism, 
the heat capacity in the K-formalism also leads to the classical gravothermal catastrophe 
(Cy < 0) in the range 0 < k < 2/3. It is worth mentioning that we have used the relation 
K = g — 1 which will be discussed in section uni 


B. Gravothermal catastrophe 

With this statistical mechanics motivation in mind, let us consider a self-gravitating fluid 
sphere at constant total energy E and volume V. In the case of an instability we have a 
gravothermal catastrophe and the total energy E of the fluid is 
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FIG. 3: The normalized heat capacities for Tsallis (dashed line) and Kaniadakis formalisms (solid 
line) as functions of k {k = q — 1). 


E = ^NkeT + p(f)$(f)d\, (5.9) 

where the second term represents the gravitational contribntions. In the light of the resnlts 
obtained in Section |Vl we can write Eq. fl5.9p as 

E = ^Nf{K)kBT +^J p{fMf)d\ (5.10) 

and nsing Kaniadakis eqnipartition theorem we can make G —> f{K)Gk in fl5.10p we have 
the gravitational potential becomes 

$(f) = -Gk [ . (5.11) 

Snbstitnting this potential into Eq. (Ih.lOp we have that the total energy of the flnid in the 
Kaniadakis formalism is given by 


E 


^Nf{K)kBT 


1 G 

27M 


Pjr) p{E) 


d^rd\' . 


(5.12) 


where we can clearly that the k- formalism fnnction /(k) affect both terms in this eqnations 
in a way that we cannot write something like /(k) 

From fl5.12p we can End an expression for the temperatnre as a fnnction of the density 
distribntion at a constant total energy. So, in Eq. fl5.12p . if E = constant, T = Tq + Ti 
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and p = Pn + pi, where the indice 0 means “nnpertnrbed,” we have the expression for Ti as 
being |35[ | 


Ti = - 


/ pi(r)<ho(r)(i^r 


(5.13) 


lNf{K)kB 

where d>o is the gravitational potential of the nnpertnrbed density distribntion (po)- Con- 
seqnently, we can directly see that the pertnrbative temperature is a consequence of KS 
analysis. 


VI. RELATIVE DEVIATIONS OF THE MODIFIED GRAVITATIONAL CON¬ 
STANTS 

In order to study the behavior of the modified gravitational constants, Eqs. fl3.2p and 
dH, when we change the k parameter, we define relative deviations functions as 

G — Gq 3 

6Gq = —-— = -K, (6.1) 

SG^ = ( 6 . 2 ) 

where we have used Eqs. fl3.2p and (13.51) in (16.11) . In Figure 4 we plot the relative deviations, 
Eqs. (16.11) and (16.21) . as functions of the deformation parameter k. Then, from Figure 4 we 
can observe that for the same k parameter the value of (16.2p is greater than (16.1|) except for 
K = 0 and K = 2/3. 

VII. CONCLUSIONS AND PERSPECTIVES 

In this work we have derived the equipartition theorem for the Kaniadakis statistics. Our 
procedure is based on the generalized Maxwellian formulation for ^-statistics. In the hrst 
application, we have analyzed the self-gravitating system in the context of Kaniadakis’ for¬ 
malism, and likewise, in TT formalism we have also obtained the gravothermal catastrophe 
for all K in the range 0 < k < |. 

In the second one, we have discussed the Kaniadakis equipartition law of energy in the 
framework of Verlinde’s entropic formalism. We have derived a gravitational constant as 
a function of k parameter. Comparing with the gravitational constant obtained by TT 
formalism, we have verified that the difference between both gravitational constants is not 
significant (figure 3). This result indicates that, in the context of the entropic framework, 
Tsallis’ formalism is more suitable when compared with k statistics because the algebraic 
expression of equipartition theorem, and consequently the modified gravitational constant 
obtained by Tsallis’ formalism is more simple. Last, the relative deviations of the modihed 
gravitational constants were studied. 

As a perspective we can analyze the formalism developed here, i.e., the connection and 
the extension of the statistical frameworks, as well as the properties of the modihed G, to 
treat high scales cosmological structures. It is an ongoing research and will be published 
elsewhere. 
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FIG. 4: The relative deviations of the modihed gravitational constants obtained by Tsallis(dashed 
line), Ea. (|6.1l) . and Kaniadakis statistics (solid line), Ea. (l6.2l) . as functions of k. 
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